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Abstract 

Within the context of Supersymmetric Quantum Mechanics and its related hierarchies of in- 
tegrable quantum Hamiltonians and potentials, a general programme is outlined and applied 
to its first two simplest illustrations. Going beyond the usual restriction of shape invariance 
for intertwined potentials, it is suggested to require a similar relation for Hamiltonians in 
the hierarchy separated by an arbitrary number of levels, N. By requiring further that these 
two Hamiltonians be in fact identical up to an overall shift in energy, a periodic structure is 
installed in the hierarchy of quantum systems which should allow for its solution. Specific 
classes of orthogonal polynomials characteristic of such periodic hierarchies are thereby gen- 
erated, while the methods of Supersymmetric Quantum Mechanics then lead to generalised 
Rodrigues formulae and recursion relations for such polynomials. The approach also offers 
the practical prospect of quantum modelling through the engineering of quantum potentials 
from experimental energy spectra. In this paper these ideas are presented and solved explic- 
itly for the cases iV = 1 and N — 2. The latter case is related to the generalised Laguerre 
polynomials, for which indeed new results are thereby obtained. At the same time new classes 
of integrable quantum potentials which generalise that of the harmonic oscillator and which 
are characterised by two arbitrary energy gaps are identified, for which a complete solution is 
achieved algebraically. 
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1 Introduction 



Factorisation methods to solve the Schrodinger equation are almost as old as quantum mechanics 
itself. Nowadays this approach is an integral part to the techniques of Supersymmetric Quantum 
Mechanics (for reviews and references to the earlier literature, see for instance Refs.fT, '2]). Even 
if only for systems with a single degree of freedom, this is a field which still offers most tantalising 
and fascinating perspectives indeed, as witnessed by a renewed strong research activity over the 
last decade at least[3], which goes on unabated and keeps exploring quite many open avenues, 
while extensions to more degrees of freedom, including fermionic or spin ones, is a realm remaining 
largely unexplored until now (for a review and further references, see Ref.[4J). 

For the single degree of freedom case which is the topic of the present work, the situa- 
tion may be characterized as follows (assuming that the potential energy is such that the energy 
eigenspectrum be bounded below but unbounded above while also being discrete and countable 
infinite without any degeneracy). Through the factorisation of the second order differential sta- 
tionary Schrodinger equation into the composition of two first order differential operators which 
are adjoints of one another, a given Hamiltonian whose energy eigenspectrum of states would be 
known, is seen to belong to an infinite hierarchy of successive intertwined pairs of Hamiltonians 
of which the energy eigenspectra may readily be identified starting from that of the first Hamilto- 
nian. All of these Hamiltonians share an identical infinite spectrum of energy eigenvalues except 
for the lowest lying state which is removed as one moves from one member in the hierarchy to 
the next. In other words, given a single degree of freedom quantum system of which the energy 
eigenspectrum is known explicitly — namely both the values of its (discrete) energy eigenvalues 
and the corresponding quantum states, the latter say in terms of their (configuration space) wave 
functions — , there arises a semi-infinite hierarchy of integrable single degree of freedom quan- 
tum systems all of whose energy eigenspectra are known likewise. Supersymmetric Quantum 
Mechanics and factorisation of the Schrodinger operator thus provide an important insight lead- 
ing towards the construction of nontrivial integrable quantum potentials, starting from already 
known ones. Given that energy eigenstate wave functions possess a number of nodes equal to the 
order of the energy level while also defining an orthonormalised countable basis of Hilbert space, 
the identification and classification of such integrable hierarchies of quantum systems should also 
provide generally for new insights into properties of orthogonal polynomials. 

However with the exception of a few cases, integrable quantum Hamiltonians are difficult 
to identify, to be used then as the starting member of such integrable hierarchies. Further re- 
strictions are required to enable the actual construction of such hierarchies. For instance if two 
successive intertwined members of the hierarchy are somehow related to one another, the induced 
recursion relations should allow for a solution which then identifies an integrable hierarchy of 
Hamiltonians. This is the basic idea of "shape invariance" [H [2] . If the potential energies associ- 
ated to two successive Hamiltonians in the hierarchy are related simply through a redefinition of 
their defining parameters, such a property of "shape invariance of the potential energy" allows for 
an explicit resolution of the energy eigenspectra of two such intertwined Hamiltonians or potential 
energies, hence the identification of an integrable hierarchy of quantum Hamiltonians. Specific 
classifications of shape invariant potentials have been achieved in the literature[Tl [2], while the 
known list — which includes further cases discovered since when the latter two publications have 
appeared^ — is certainly far from exhausting a complete classification (for recent developments 
and further references to the relevant literature, see for instance Refs.[5l El [TJ El [9j). 

The motivation and results of the present work follow a similar logic, but rather than 
compare two successive intertwined Hamiltonians in the hierarchy, it considers what may occur 
when a Hamiltonian, or potential energy, further up in the hierarchy than simply the next one, is 
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related to the first one in the hierarchy. More specifically, we shall be considering the restrictions 
arising for the hierarchy and the potential energies defined by it, when the ground state wave 
functions of two members of the hierarchy happen to coincide. This is a condition stronger than 
requiring shape invariance for the corresponding two potential energies. Indeed, as is well known, 
from the ground state wave function one readily identifies the potential energy, and thus implicitly 
also the entire energy spectrum. Hence in actual fact, if two Hamiltonians in the hierarchy share 
a common quantum ground state, up to a constant shift upwards in their energy spectra they 
share an identical potential energy and thus energy eigenspectrum inclusive of the wave functions 
for all their energy eigenstates. As a consequence, and in a manner similar to what happens when 
imposing shape invariance of the potential energies, it becomes possible to construct an integrable 
hierarchy of quantum systems, leading further to specific properties for the associated orthogonal 
polynomials. Clearly the above restriction may somewhat be relaxed by requiring only shape 
invariance rather than an exact identity for two potential energies defining two Hamiltonians 
separated by > 1 steps or levels in the hierarchy, the case = 1 being the only one heretofore 
considered in the literature in the context of shape invariant potentials. 

Leaving aside for future work the general case, the present paper details the consequences of 
having a common quantum ground state for a given Hamiltonian and the next-to-next one in the 
hierarchy, namely the case N = 2. This is done by specifying the energy gaps between the first 
two energy levels of the starting Hamiltonian, with as a result the energy eigenvalues simply being 
the repeated duplication of these two gaps thus producing a two gap periodic energy spectrum. 
Clearly the harmonic oscillator is one particular degenerate case of this general situation when 
the two energy gaps are identical, with as consequence that two intertwined Hamiltonians in 
the hierarchy then share the same quantum ground state while the spectrum simply becomes 
equally spaced. The latter simple case will serve the purpose of a warm-up illustration for the 
more general analysis which follows when requiring identical quantum ground states for two 
Hamiltonians separated by A^ levels in the hierarchy. The cases A^ = 1 and N = 2 will explicitly 
be solved in this paper. The generic situation for A^ > 3 is left to future work which is still 
on-going. 

Obviously, the general case for arbitrary A^ opens the way towards quantum potential 
engineering to produce quantum systems with A^ gap periodic energy spectra of which the A^ first 
energy gaps take specific values, predetermined for instance from experiment. The approximation 
of actual physical systems provided by such engineered quantum models is expected to grow 
better and better as the value of A^ keeps increasing, in a manner dependent of the values of 
the A^ first physical gaps to be reproduced. Furthermore, it is worth noting as well that such 
quantum potential engineering must also be related [3j to the well established field of inverse 
scattering methods in integrable Hamiltonian systems and hierarchies which allow in principle 
for a reconstruction of a quantum potential from the knowledge of its energy eigenvalues and 
phase shifts (see for instance Ref.[10j). Finally, given the periodic structure of such hierarchies 
involving a finite number of characteristic classes of wave functions and associated orthogonal 
polynomials, specific recursion relations between classes of orthogonal polynomials ought to be 
generated from such constructions. In these two latter contexts as well, one could thus also expect 
new developments following from Supersymmetric Quantum Mechanics. 

The paper is organized as follows. Section 2 briefly recalls the main features of Super- 
symmetric Quantum Mechanics for a single degree of freedom system on the real line, and how 
hierarchies of integrable quantum Hamiltonians may be constructed out of a Hamiltonian whose 
energy eigenspectrum is known. Based on this understanding the programme described above 
in terms of periodic hierarchies is outlined. Section 3 then applies these ideas to the simplest 
situation, with A^ = 1, to show that the ordinary harmonic oscillator is recovered while known 
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properties of the associated Hermite polynomials are reproduced using the methods provided 
by Supersymmetric Quantum Mechanics. Then in Section 4 we turn to the actual original con- 
tent of the present work, by solving completely the N = 2 periodic hierarchy. Given the two 
gaps characteristic of that case, the corresponding potential energies are determined, as well as 
the energy eigenvalues and eigenwave functions for all states in the hierarchy. It is established 
that the associated ensemble of orthogonal polynomials are the generalised Laguerre polynomials. 
Two new generalised Rodrigues formulae are thereby discovered for these polynomials, while the 
methods of Supersymmetric Quantum Mechanics provide for specific recursion relations for these 
functions. Some interesting properties related to the singularities of the corresponding potentials 
are also commented on. Finally in Section 5 we present our main conclusions. 



2 Engineering Hierarchies of Integrable Quantum Hamiltonians 

2.1 Some Basics of Supersymmetric Quantum Mechanics 

Also for the purpose of establishing our notations and conventions, let us briefly outline the 
salient features of a given hierarchy of intertwined quantum Hamiltonians, without providing 
all the justifications — readily available from the literature[Tl [2| — for the results stated hereafter. 
Levels in such a hierarchy are labelled by an index A = 1,2, ••• corresponding to quantum 
Hamiltonians Hx, Hi for A = 1 being the first member of the hierarchy. In the configuration 
space representation of quantum states in terms of real wave functions ip{x) with, in the present 
work, the configuration space variable x assumed to take values in the real line, x G M, the 
quantum Hamiltonian operators are expressed as, 

Vx{x) being the corresponding potential energies for the single degree of freedom x of mass mg. 
The stationary Schodinger equation determines the energy eigenspectrum of each system in the 
form, 

Hx-^xM = Ex,n ^\,n{x). n = 0, 1, 2, • • • . (2) 
The eigenwave functions are assumed to be orthonormalised, 

dx %Ij*x^^{x) %l}x,m{x) = 5n,m, 71, m = 0, 1, 2, • • • , (3) 

) 

while for later purposes it is useful to consider the successive gaps in each of these energy eigen- 
spectra, 

AA,n = ^A,n - ^A,n-l > 0, n = l,2,3,---. (4) 

For convenience, let us also introduce the following change of configuration space variable, u G M, 
such that, 



the variable u thus having the physical dimension of the inverse square root of energy, [u] = E^^/"^ . 
The above relations then take the form, 
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(with a slight customary abuse of notation for these mathematical functions). 

Factorisation of the stationary Schrodinger equation is achieved as follows. The ground 
state wave function of each Hamiltonian uniquely determines a superpotential W\{u) through, 

Note that this relation integrates to, 

V'A,oH = iVA,oe-i''^"^^("), (8) 

Nxfi being a normalisation factor (which may be chosen to be real and positive). In terms of 
Wx{u), one then constructs two linear first order differential operators which are adjoint of one 
another and factorise the Schrodinger equation. 



d _ . d 

and such that. 



H,=A{Ax + E,,„, V,(u) = Wi{u) - W{(u) + (10) 

where W'^^u) = dW\{u)/du (hereafter the notation du = d/du is often used as well). 

Two intertwined Hamiltonians H\ and -ffA+i belonging to a same hierarchy are then such 

that, 

^^A+l = ^1+1 ^A+l + ^A+1,0 = ^A ^1 + ^A,0, (11) 
with as consequence the following intertwining relation for the corresponding superpotentials, 

Vx+i{u) = Wl^, - W'x+, + Ex+i,o = Wl + W'x + Ex,o. (12) 

Generally there is no known method for solving any of these two types of Riccati equations, 

- + Ex,o = Vx, - W'x+i + Ex+1,0 = + W'x + Ex,o, (13) 

for the superpotentials Wx{u) and Wx+i{u) given, say, the functions Vx{u) and Wx{u), respec- 
tively. However knowledge of the ground state wave functions V'A,o(^) and V'A+i,o(^) provides 
just such a solution. As a matter of fact, knowledge of the whole energy eigenspectrum of Hx 
provides at once not only for the superpotentials Wx{u) but also for the whole energy eigenspec- 
trum of -ffA-i-i) thus in particular also Wx+i{u) from the ground state ipx+i,o{u)- Indeed, the 
following relations apply between the energy eigenspectra of the two intertwined systems which 
are mapped into one another by the operators Ax and except for the ground state of Hx which 
is annnihilated by Ax (hereafter, n = 0, 1, 2, • • • ), 

tpX+l,n{'^) = -y=^=== ^A V'A,n+l(«), i^X,n+liu) = == A\ V'A+l,n(«), (14) 

V-^A,n+l — ^X,0 \/-l^X,n+l — ^X,0 

Ex+l,n = Ex,n+1, (15) 

while, by construction, 

Ax^x,o{'^) = ^- (16) 

Consequently any such semi-infinite hierarchy of pairwise intertwined Hamiltonians is char- 
acterised by two complementary sets of information. On the one hand, an ensemble of energy 
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gap^ A;^=i „ = hui^n = 1;2, • • • ), and on the other hand, an ensemble of intertwined super- 
potentials, Wx{u) (A = 1, 2, • • • ), obeying the recursion relation^, 

Wl^,{u) -W;^^,{u) + Ax^i = Wl{u) + W'^{u), Aa,i=Ai,a, A = 1,2,.... (17) 

Any eigenstate ipx,niu) (n = 1,2,---) of Hx may then be constructed from the ground state 
ipx+nfii"^) of Hx+n through the repeated application of operators defined in terms of the 
superpotentials, 

= [{Ax,n + Aa,„_i • • • + Aa,i) • (Aa,„ + Aa,„_i + • • • + Aa,2) • • • (AA,n + Aa,„_i) • Aa,„]"^/2 X 

X 44+l---4+n-lV'A+n,oN, n = 1, 2, • • • , (18) 

while all ground state functions are themselves directly constructed from the superpotentials, 

V'A,o(n) = iVA,oe-^^"^^("\ (19) 
with the normalisation factors Nxfi remaining to be determined. 

2.2 Engineering of Quantum Energy Spectra 

It may so happen that two members of the hierarchy, say Hx and Hx+n with > 1, share a 
common ground state, or equivalently a same superpotential, Wx+Ni^) = Wx{u). Since these two 
systems then share the same potential energy up to an overall shift in the ground state energy of 
Hx+N given by the sum of the first N gaps in the spectrum of Hx, a periodic structure of order 
N exists in the hierarchy which then also readily extends to — oo < A < +00. 

Consequently without loss of generality one may consider that the Hamiltonians which 
share the same ground state are Hi and i^i+Ar. Such an infinite hierarchy is then constructed 
(up to the arbitrary energy level -^1,0) from the knowledge of A'^ gaps, namely those o^ Hi, 
An = Ai^n = fi^n for n = 1, 2, • • • , A', and A^ superpotentials, namely Wx{u) with A = 1, 2, . . . , A^ 
since Wi+n{u) = Wi{u), while all these finite number of data must solve the A^ recursion relations 
which close back onto themselves since we now also require that Wi+^iu) = ^^1(^)5 

W^^iiu) - W'^_,i{u) + Ax = Wliu) + W'^iu), A = 1,2, • • • ,iV. (20) 

Because of the isospectral properties of the hierarchy except for its ground states which 
are removed or added as one moves up or down in the hierarchy using the operators Ax and 
A^, an order A^ periodicity arises. Given any Hamiltonian Hxq, all its partners Hx at the levels 
Aq (mod A^), namely with A = Aq + kN (/c G Z), share the same energy spectrum up to a shift in 
energy by A;(Ai + • • • + Ajy), and share the same energy eigenwave functions, V'A,n(^i) = V'Ao,n(^) 
(n = 0, 1, 2, • • • ). Furthermore, except for their ground state energies, the gaps of each of these 
spectra consist of the periodic duplication of their first A^ gaps while preserving their order. 

^Note that since AA+i,n = ^\,n+i {n = 1,2, ■•■), it suffices to known the energy gaps Ai,„ of the first 
Hamiltonian of the hierarchy, Hi . All energy spectra are known from this spectrum of gaps up to the arbitrary 
overall shift by ground state energy of the first hierarchy, -Ei.o- 

^Note that from this point of view, provided these Riccati equations may be solved, it even becomes possible 
in principle to extend the hierarchy to negative values of A < 0, hence extend the semi-infinite hierarchy to a truly 
infinite hierarchy of integrable quantum Hamiltonians. Applying this remark to the harmonic oscillator as the 
choice for Hi has produced some interesting new infinite classes of integrable quantum svstems[5llllj. 

^Note that Aa.^ = Ai,„+a-i. 
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Finally, as one considers each of the Hamiltonians Hi to Hj^ in turn, their first gaps consist 
of the N cyclic permutations of the N gaps (Ai, A2, • • • , Aat) in that order, with Ai being the 
first gap in the spectrum of Hi, A2 in that of H2, and so on, until H^^i of which the first gap 
is again Ai. 

Such a specific situation thus possibly opens the avenue towards engineering quantum po- 
tentials of which the first energy levels take prespecified energy values, for instance dictated 
from experiment. The corresponding energy spectra are then gap periodic. Such models may 
provide good approximations to actual physical quantum systems, with a quality of approxima- 
tion that presumably would improve as the value of increases. Once the first A^ gaps are 
specified, there remains to solve the A^ coupled nonlinear first order differential Riccati recursion 
relations ()20p . for which even a numerical approach could be developed for practical applications. 
Work on the generic case with A^ > 3 is still on-going. Hereafter the solutions for the cases A^ = 1 
and N = 2 are detailed. In a certain sense the equations to be solved then define a degeneracy 
of the general case when A^ > 3. 



3 The TV = 1 Case 

The simplest case A^ = 1 is characterised by a single gap and a single superpotential. Up to the 
value of their ground state energy which is shifted by a multiple of that gap, all Hamiltonians 
of this hierarchy, of periodicity N = 1, thus define the same quantum system with an equally 
spaced energy spectrum. Quite obviously, all these Hamiltonians are those of the same harmonic 
oscillator whose angular frequency, ojq, is set by the choice of unique energy gap. Using the 
notations 

Ai = A = ;;^o, Wi{u) = W{u), ^1,0 = ^0, (21) 
the single Riccati "recursion" relation to be solved is 

W^{u) - W'{u) + A = W^{u) + W'{u). (22) 

Clearly there always exists a solution, given as, 

W'{u) = ^A, Win) = l^{u- no) , (23) 

uq being an arbitrary integration constant. Correspondingly one finds, 

Vi{u) = -W + Eo = ^A\u-uo)^ -^A + Eo = ^mooj^oix-xo)^ -^hwo + Eo, 

V2{u) = W^ + W' + Eo = ^A\u-uof + ^A + Eo = ^moio^o{x-xof + ^nu;o + Eo,{2A) 

which are indeed two choices of potential energies for the harmonic oscillator of angular frequency 
Wo and centered at x = xq, shifted by its quantum energy gap A = hujQ. 

Given p^ . the ground state wave function inclusive of its normalisation is readily found 

to be, 

V'io(n) = (^y^" e-i^(— = (H^y^' e-^(— = (H^y^' e-^^^\ (25) 
where the following change of (dimensionless) variable is introduced, 



^[n-uo) = J'^{x-xo)eR. (26) 
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The energy spectrum of the first Hamiltonian in the hierarchy, Hi, is simply, 

Ei^^ = Eo + nA, n = 0,l,2,---, (27) 
while the energy eigenwave functions are constructed as, given (|18p . 

V'i,n(^i) = ^=1== (^t)" V'i,o(u), n = 0,l,2,---, (28) 
V A" • n! ^ ^ 

where one has, 

A = Ai = du + W{u) = ^{v + d,), = A\ = -du + W{u) = ^{v - d,). (29) 

Consequently, one readily finds, 

i^i,n{^) = i^T^) -j==e-^^ Pn{v), n = 0,l,2,---, (30) 

where Pn{v) are polynomials in v of order n defined by the formula, 

Pn{v) = {v - dX e--^''", n = 0,l,2,..., (31) 
and obeying the following orthonormality properties, given the conditions in ([6]), 

/ dve-'' Pn{v)P„^{v) = 5n,m2''n\^/^, n,m = 0,1, 2,--- . (32) 

J —oo 

That these polynomials are precisely the usual Hermite polynomials should be quite obvious. For 
instance using the fact that 

e-i-^ (i;-4) = (-4)e-^', (33) 
implies that the above generating formula for the Pn{v) polynomials is reduced to, 

P„(7;) = e^%-^' {v-d,T e-^"' = {-If e^" e"-' = i/„(t;), (34) 

which is indeed Rodrigues' formula defining the Hermite polynomials, which also obey ([32]) (see 
for instance Ref . 1121 ) . 



This identification may also be achieved by considering the original Schrodinger equation, 
say for Hi , in which the known potential energy, energy spectrum and wave functions are substi- 
tuted, namely, 

-dl + iA2(n - uof - + ^i,„(n) = (Sq + nA) V'i,n(n), (35) 
or equivalently, 

{dl - v^) ^i,„ = -(2n + 1) Vi,n. (36) 
In terms of the polynomials Pniv), one then finds, 

d 



dv^ 2.- + 2„jP„(„)=0. (37) 

This is indeed the second order linear differential equation of which the normalisable solution is 
proportional to the Hermite polynomial Hn{v) (the other linearly independent solution not being 
normalisable over M)|12j. 
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Given the fact that the operators A and map these different quantum states into one 
another, specific recursion relations may also be established for these polynomials by exploiting 
the super symmetric quantum mechanics realisation of the hierarchy. For instance the relation 
(jl4p which in the present situation becomes, 

tpl,n+l = ^==== 'lp2,n, '4'2,n = tpl,n, n = 0, 1,2,---, (38) 

^/{n+ 1)A 
translates into the relation, 

Prr+liv) = ei-' {V - 4) (e-^' Prriv)) = 2vPn{v) - (39) 

However by taking the derivative of this identity and with the use of the differential equation 
(j37|) . we also have, 

- 2vPn^l{v' 



dv dv \ dv 

= 2nP„_i(t)), (40) 
SO that the above recursion relation also writes as, 

Pn+l{v) = 2vPn{v) - 2nPn-l{v). (41) 

Given the initial values Pq{v) = 1 and, by definition, P-i{v) = 0, this is indeed the three steps 
defining recursion relation for the Hermite polynomials [12], displaying once again explicitly how 
the methods of Super symmetric Quantum Mechanics provide insight also into the properties of 
orthogonal polynomials. In the present case of course, nothing new in this respect is gained from 
the present approach for as simple a system. 



4 The N = 2 Case 

The hierarchy corresponding to the case N = 2 displays a periodicity of order N = 2. It is 
characterised by two independent gaps, Ai = huji and A2 = huj2, and two superpotentials, Wi{u) 
and W2{u), with the restriction that W3{u) = Wi{u) when considering the Riccati recursion 
equations, 

Wi-W2 + Ai = Wl + W[, 

Wl-W[^^2 = W^ + W^. (42) 
Taking the sum and the difference of these two equations, one finds, 

W^ + W{ = ^{A2 + Ai), Wi -Wl = ^{A2 + A^), (43) 

of which the general solution is obviously, 

W,iu) = hA2 + Ai){u-uo)-l^'~^' ^ 



4' 2A2 + Ai u-uo' 

W2{u) = \{A2 + Ai)(n - ^o) + ^ ~ (44) 

4 Z L\2 + iAi U — Uq 
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Correspondingly the two potential energies of the hierarchy are obtained as, 

. l(A..A,^(.-„„)^.l '^'-^^;f^^y ^J^^^ ,}.^,,,,,. (45) 

Note how as is indeed expected, the two quantum systems Hi and H2 are transformed into one 
another under the permutation Ai •(-)• A2 and an upward shift in energy by Ai for the energy 
spectrum of H2 as compared to that of Hi . Another remark which provides a useful check on the 
results hereafter, is that the particular degenerate choice A2 = Ai indeed reproduces exactly the 
previous case = 1 as it should, namely the ordinary harmonic oscillator. 

These two potentials are even in (u — uq). Unless we have the degenerate situation with 
A2 = Ai, both potentials are singular infinite at = uq, one running to (—00) and the other 
to (+00) in a symmetric fashion around u = uq, depending on the sign of (A2 — Ai). Hence 
one of these potentials is then not bounded below, and yet this bottomless throat-like potential 
well remains sufficiently narrow so that the spectrum of energy eigenstates is bounded below 
and infinite discrete with normalisable wave functions in that case as well, as the construction 
hereafter establishes explicitly. 

4.1 Energy eigenspectra 

For notational convenience it is useful to introduce the following parameter measuring the relative 
difference of the two gaps, 

1 A2-A1 1^1 1 A2 1 Ai 

a = — — , < a < -, — ha = - — , a = • 4d 

2A2 + A1' 2 2' 2 A2 + A1' 2 A2 + A1 ^ ^ 

The value a = then corresponds to the = 1 hierarchy of the harmonic oscillator, while the 
two quantum systems Hi and H2 are simply transformed into one another by changing the sign 
of a. From the two superpotentials 

Wi{u) = hA2 + Ai){u-uo) W2{u) = ]{A2 + Ai){u-uo) + (47) 

4 U — Uq 4 U — Uq 



and the construction (jl9p . one readily finds for the two ground states, 
^i_o(u) = Afi,o|n-no|"e-i(^2+Ai)(u-uo)2^ ^^^^^^^ ^ ^^^^^^_^^^-a ^-l(A^+Ai)iu~uo)^ ^ (48) 

with their normalisation factors given by. 



1/4 

A^i^o = A^o , ^ N2,o = No , ^ No = ( "^"^ ^ ) . (49) 



1 .r _.r ^. _ f mo{A2 + Ai) 



ra + a) A/r(i-a) 



2/12 

This form of the solutions for the ground state wave functions invites the following change 
of (dimensionless) variable, 

. = ^(A2 + Ai)iu - no? = \ H^^i^llM (, _ ,„)2. (50) 

However care must be exercised when using the inverse relations for this change of variable, since 
two separate situations — to be distinguised by the notations v+ and v- when necessary — must 
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then be considered, depending on the sign of {u — uq), namely, 



u-uo>0 : u-uo = H — , ^/v^, v+ > 0, 

V^2 + ^1 

2 

u — uq <0 : u — uq = = ^/vZ, ^;_ > 0. (51) 

V A2 + Ai 

Most if not all expressions listed hereafter apply only in the domain (u — uq) > and for the 
variable f+, even though this fact will then not be emphasized explicitly and only the notation v 
be used. When extending expressions into the (n — uq) < domain, one must beware of possible 
changes of sign induced by the above change of variable in terms of v-. Given these cautionary 
remarks, the ground state wave functions are given as (these expressions apply as such whether 
in terms of v+ or V-), 

Mu) = No r-i/2(i/2 + a) e--2\ V2,o(^^) = Nq T-^/^{1/2 - a) 6"^^ (52) 



Through the construction in (jlSp . all other energy eigenwave functions may readily be 
constructed through the repeated application of the A\ and A2 operators on these two ground 
state wave functions. Furthermore, given their two gap periodicity the energy spectra of both 
Hamiltonians Hi and H2 are easily identified. One finds, with p = 0, 1, 2, • • • , 

Hi : Si,2p = (A2 + Ai)p + Ei,o, 

^i,2p+i = (A2 + Ai)p + Ai+Ei,o = (A2 + Ai)(p + ^-a)+Ei,o, (53) 

H2 : S2,2p = (Ai + A2)p + Ai + Ei,o = (A2 + Ai)(p + ^-a)+^i,o, 

£^2,2p+l = (A2 + Ai)(p + 1) + ^1,0. 

Given the above change of variable, one now has, 

A{ = Ta^TaT - ia^-V2 _ „i/2^^^ = Ta^TaT + ^-av'^/^ - v'/'d,^ , 

(54) 

as well as, 

A, = Ta^TaT Qz;1/2 - ia^;-V2 + ^V^^,^ , = v^A^TaT Q^V2 + 1 «^-i/2 + ^1/2^^ 

(55) 

Note well that these expressions apply in the u+ domain. In the V- domain, an extra minus sign 
multiplies each of the expressions in the r.h.s. side of these four identities. Consequently all energy 
eigenwave functions obtained from the action of an even number of A\ and A2 operators possess an 
identical expression whether in the f+ or the V- domain, while all those wave functions obtained 
from the action of an odd number of these operators possess expressions of opposite sign in these 
two domains. This observation plays an important role in the orthonormality properties of these 
quantum states. Note also that one may expect that the repeated application of these operators on 
the wave functions defining the two ground states would produce terms with powers of v growing 
more and more negative as the number of operators increases, given the last two contributions, in 
u"^/^ and v^^'^dy, in the above differential operators. However as will be established, cancellations 
are such that this never happens, and indeed as ought to be expected, up to some overall factor 
function of v including the exponential factor e"""/^ and a power of v, the wave functions are 
given by some polynomial in v of finite order. 
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To see how this comes about, let us first compute some of the lowest lying states, beginning 
with the Hi system. A direct calculation finds, 

= ^4V2,oH = iVo^-l/2(3/2-a)^;-^"+5e-^^ 
v^i 



L== 44 V'l.oW = iVo r-i/2(3/2 + a) e-"2- (v-\-o), 



V'i,3M = , , ^ 444v^2,oW (56) 

VAi(Ai +A2)(Ai + A2 + Ai) 



Likewise for i?'2, a direct calculation of ^'2,1(^)5 ^2,2(1*) and ■02,3 (i*) using the same approach 
finds the same final expressions as above with simply the substitution of a by (—a) everywhere, 
namely, 

0^2,1^ = ^4Vi,o(n) = Ao^-V2(3/2 + a)^;l"+le-|^ 
V A2 



02,2(n) = ^=^==AlA\^2fi{u) = NoT'^l\-i/2-a)v~'2'-e-'2- (v-\ + a 
vAi(Ai + A2) V ^ 



^^2,3^ = ^ , , :444V^1,0(U) (57) 



^A2(A2 + Ai)(A2 + Ai + A2) 



In these calculations the cancellation mechanism mentioned above is already operational starting 
with the second excited level, n = 2. This is a general feature which extends to all excited states. 

A careful analysis for the other energy eigenstates uncovers the following structur^. For 
the Hamiltonian Hi, one finds, for p = 0, 1, 2, • • • , 



^/p\T{p + 1/2 + a) 



0i,2p(n) = , =z;^°e-2-A7°(^), 



V'i,2p+i(n) = , . ^ ^^"e-^A°f(.), (58) 

A/p!r(p + 3/2 - a) 



where the polynomial factors are obtained from, 

P{=™^(i;) = ^;-|°e^''^)P(a)^;^"e-^^ 



_i /o 1 /o , \ — ia —if 



1 1 „ 1 . 



= ^;-2 + 2"e2^L>P(a)^;2-2"e-2^ (59) 
in terms of the nonlinear second order differential operator D{a) defined by, 

D{a) = Qt;V2 _ 1^^-1/2 _ ,1/2^^^ Q^i/2 ^ 1^^-1/2 _ ^1/2^^^ _ (go) 



^Once again let us recall that these expressions apply in the «+ domain. In the V- domain, those for i/)i,2p+i(ti) 
and V'2,2p+i('") acquire an overall additional minus sign, while those for i/'i,2p(^t) and i/'2,2p(w) remain unaffected. 
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Likewise, a direct analysis confirms again the {a -f-)- —a) substitution rule for passing from the 
Hi system to the H2 one. One finds, for p = 0, 1, 2, • • • , 

Vp!r(p+ 1/2 - a) 

V'2,2p+i(n) = ^==^^===v\^\'^ e-\^ P^2f(v\ (61) 
A/p!r(p + 3/2 + a) 

where, 

= v-^2-¥e2" DPi-ojv^^+l'^e-^^^ (62) 
The previous explicit examples thus correspond to, 

Pf,r = l, Pll--{v)=v-\-a ■ P2T(^;) = 1, p-y{v) = v-\ + a, 

P'^f{v) = l, P^f{v)=v-\ + a ■ P^f{v) = l, p-f{v) = v-l-a. (63) 

Let us compare these expressions with those for the generalised Laguerre polynomials defined by 
their Rodrigues formula, 

= - x-^ — e--) , n = 0, 1, 2, • • • , (64) 



with in particular, 

L^^\x) = l, 4^)(:r) = (l + 7)-x, (x) = ^(7 + 2)(7 + 1) - (7 + 2)x + Jx^. (65) 

Hence we may certainly write. 



P?f{v) = P?f{v) = -L^r'^\v), (66) 

and likewise for P|,o'"'(^)' P^^iv). P2fl^{v) and Pl^f'iv) under the substitution (a ^ -a). These 
observations thus suggest the following general identificatiorH (in the vj^ domain), 

PtT{v) = i-irpl4-'^^"\v), P°dd(^) ^ i-l)Pp\Ll^-''\v), (67) 

ATW = (-l)^p!4"^""^^), ^^2,^ (^) = {-irp\Ll^'"'\v). (68) 

That these identifications are indeed correct will be confirmed hereafter. Assuming this to 
be the case, in conclusion we have thus obtained for both Hi and H2 their spectra, see (f53l) . but 
also their energy eigenwave functions, in the form. 



^The factor p! and the sign ( — 1)'' follow from considering the factor multiplying the highest power in v as 
generated for both classes of polynomials from their respective formulae above. 
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where A'^q = (mo(A2 + Ai)/(2/i^)) In particular, note how in the hmit when A2 = Ai, namely 
a = 0, all these results do reproduce those for the harmonic oscillator of angular frequency 
ujQ = (A2 + Ai)/(2/i), given the identities [T^], 

H2p{x) = (-If 22>!4~^)(^2)^ i/2p+i(x) = (-1)^22^+1 p!x4^^(x2). (71) 
That the orthonormality properties ([6]) are also obtained follows from the identity [12]. 

t-OO 1 

/ dxe-'''x^L'^'Hx)L'^\x) = 6n,m-,nj + n + l), n, m = 0, 1, 2, • • • . (72) 
Jo nl 

For instance for the Hi system (and likewise for the H2 system), through a careful analysis of the 
contributions from the V-^- and V- domains the normalisation of the states V'i,2p(i*) and V'i,2p+i(^i) 
requires, respectively, 

dve-''v-"^+''{l}i^'^"'\v)^ = ^r(p+l/2 + a), 

dv e-" v'—'^ {l}y"\v)^ = ^r(p + 3/2-a), (73) 

which is indeed correct. Similarly the overlaps between two "even" states V'i,2p('u) (or ip2,2piu)) on 
the one hand, or two "odd" states 'ipi,2p+i{u) (or '4>2,2p+iiu)) on the other hand, for two different 
values of p, are directly seen to be vanishing as they should, given the result ([72]) . Finally, the 
overlap between an "even" and an "odd" state for a same Hamiltonian, Hi or H2, is found to 
vanish as well, on account of the specific changes of sign that arise when considering the overlap 
integral over u G M decomposed in terms of the two domains v+ and V- . 

4.2 Orthogonal polynomials generated by the N = 2 hierarchy 



In order to determine which are the orthogonal polynomials Pl''p^{v) and Pf^ (?;) (and P2'p"[v 



and P2p'^{v)) generated by this N = 2 construction of super symmetric quantum Hamiltonian 
hierarchies, let us turn to the stationary Schrodinger equation defined by Hi (or H2), since both 
the energy eigenvalues and eigenwave functions are known. Given the potential energy Vi{u) 
determined previously, and in terms of the change of variable v or v±, the stationary Schrodinger 
equation foi|j Hi becomes, 

.2 1^ , 1 la(«-l) If , A\ I r ^ {Ei,n-Eifl) 



-vd^ - + 4^ + — 2 + 2 J J = A2 + a; 

Distinguishing now the two situations when n = 2p or n = 2p + 1 , on account of the expressions 
in (jSSp one finds for the relevant polynomials, 



®The same detailed analysis for H2 leads to the same conclusion, and amounts simply to the substitution 
everywhere of a — j- — a. 
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v-^+(l-a-v + l]-^+p]P^'l''{v)=0. (76) 



dv"^ \ 2 J dv 

Since the generalised Laguerre polynomials Ln^ (x) are the normalisable solutions (for the measure 
dx e~^x'^) to the linear second order differential equation[12j 

d \ 

^dx^~^^^~^^'^^d^^'^) n = 0, 1,2,---, (77) 

indeed, up to a normalisation factor that has been determined previously already, the polyno- 
mials Pi^™{v) and Pip'^{v) coincide with the generalised Laguerre polynomials Lp ^"'""^(w) and 

{--a) 

Lp^ (v), respectively. The previous conjectured identification is thus established. The N = 2 
hierarchy of integrable quantum Hamiltonians as constructed following the programme outlined 
in the present work is indeed a hierarchy generating the generalised Laguerre polynomials as the 
associated set of orthogonal polynomials. 

4.3 A generalised Rodrigues formula for generalised Laguerre polynomials 

Given the identification of orthogonal polynomials having been achieved from the Hi system, we 
have thus established the following formulae for generalised Laguerre polynomial^. 



S "Vv) = LlZ^-i+i"ei^L>P(a)i;t-3"e-3^ (78) 
p\ 

By introducing a parameter 7 defined to be either 7 = — 1/2 + a (given the first identity) or 
7 = 1/2 — a (given the second identity), in either case these two formulae are reduced to the 
following two generating formulae for generalised Laguerre polynomials, 

L'^^v) = D^(±7) (79) 

where the following two nonlinear second order differential operators are introduced, 

D,{±j) = - ^ ± 7) V-'/' - v'/'d^^ + i ± 7) V-'/' - j . (80) 

Note that these same operators may also be expressed as. 

As compared to the usual Rodrigues formula for generalised Laguerre polynomials, 

LW(^,) = le^^;-7 (^A.y e-''vP+\ (82) 



the two new generalised Rodrigues formulae for generalised Laguerre polynomials in (j79p are 
noteworthy. First note that irrespective of whether one uses the operator .02(7) or the operator 
1^2 (—7)5 the same generalised Laguerre polynomial is generated from (I79p . Such a very specific 



^The formulae obtained from the H2 system are the same under the substitution (a — > —a). 
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feature can only apply provided once again a series of cancellations occur when expanding the 
action of the powers of these operators. This may indeed be confirmed by a direct evaluation of 
some of the polynomials of lowest order. 

A second noteworthy property of the formulae in (j79p is the following. In contradistinction 
to the usual Rodrigues formula (|82p which makes it obvious that the p order action of the linear 
first order differential operator d/dv generates a polynomial in v of order p, the formulae in (|82p 
may naively seem to generate some polynomial in v of order 2p because of the p order action 
of the nonlinear second order differential operators 1^2 (±7) of which the leading term is in v. 
Yet this not the case, once again owing to specific cancellations that occur such that no negative 
powers of v remain either in spite of the contributions in \/v in D2{-^^), while also keeping the 
polynomial of order p m. v and not higher. 

In these two respects the formulae in (j79p are indeed totally different from the usual Ro- 
drigues formula for generalised Laguerre polynomials. To the authors' best knowledge, the gener- 
alised Rodrigues formulae in (|79p are new in the literature. And they indeed follow directly from 
the structure of the N = 2 hierarchy of integrable quantum Hamiltonians constructed along the 
lines of the programme outlined in this work. 



4.4 Recursion relations 

As a first way of generating a recursion relation between polynomials of different order in v, 
consider the formulae (I79p in the following form, 



p-t-i (p + In 



= —v~-i~^^e^''D2{±-f)(v-i+^"'e-^''L'^\v)). (83) 

p + 1 V ^ / 

An explicit evaluation of the latter expression then finds, once again following the cancellation of 
a number of terms, 

= ^ ((^ + 1 - ^) - (7 + 1 - 2^) d.L'fHv) - vdlL(^\v)) . (84) 

However by using ([77|l to substitute for vdlL^\v) in this identity, one has, 

= ^ ((^^ + 1 + ^ - ^) + ^d.L(;\v)) , (85) 

which is an already known recursion relation for these polynomials [12]. Nevertheless, this con- 
clusion confirms once again the correct identification made for the polynomials related to this 
N = 2 hierarchy construction, inclusive of their normalisation factors and choices of sign. 

Since all the energy eigenwave functions of both Hamiltonians Hi and H2 defining the 
N = 2 hierarchy are expressed in terms of generalised Laguerre polynomials, through the action 
of the operators Ai, A2, A\ and A2 on these states it becomes possible also to derive a series of 
recursion relations for these polynomials. Possibly some of these relations could be new. This 
is to be done by considering the general relations in (jl4p and using all four possible actions 
of the operators mapping between levels within the hierarchy, whether the excitation level n 
is even or odd. A detailed and careful analysis of all these possibilities leads to the following 
four independent relations, after the introduction of the relevant parameter 7 as a function of a 
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depending on which case is being considered, 

P + 1 V2 2 V2 y J 

= ^-l-he^ Q^i/2+ 1 Q+^^^-i/2_^i/2^^^ ^i+i7e-^i(7)(^^ (86) 
L%'\v) = --±-y\-he^^ Q,i/2 + 1 Q _ ^-1/2 _ ^1/2^^ j 

Working out these expressions exphcitly, once again cancellations occur and one is left with, 

4<-'\v) = -l-(^ + ^d,)L(7)(^), 

4^+1) (.) = -d,L^;Uv), 

L^^^'\v) = {l-d,)L^;'\v), (87) 
= -^{v-7-vd.)L(;^\v). 

Using for instance the second of these relations, the other three become, 

4^+^)(t;) = L'f\v) + L^;L\'\v), (88) 
= -^{iv-^)Lp\v)+vL^;_\'\v)), (89) 

while (|85|) then reduces to another three terms recursion relation, 

= ^ ((p + 1 + 7 - ^)4"^(^) - • (90) 

By taking different linear combinations of these recursion relations, one reproduces recursion 
relations for generalised Laguerre polynomials known from the literature and standard Tables [12]. 



4.5 Further comments 



Given the above complete resolution of the N = 2 case for a periodic hierarchy of integrable 
quantum Hamiltonians using the methods of Supersymmetric Quantum Mechanics, a few more 
remarks are probably useful and relevant. 

Among the examples of shape invariant potentials and their corresponding superpotentials 
listed on pages 40 and 41 of Ref.[2], the second one is that of a superpotential of the form (in 
units such that 2mo = 1 and ^ = 1), 

W{r) = Lr -^-±1, V{r) = ]u\^ + ^-^^-(i + l)co, Eo,i = 0, (91) 
2 r 4 \ 2 / 

r > being the radial coordinate of the three dimensional spherically symmetric harmonic os- 
cillator, and £ = 0, 1, 2, • • • the usual orbital angular-momentum quantum number, leading to an 
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equally spaced spectrum of eigenvalues, Ei^n = 2nw (n = 0, 1, 2, • • • ). With the correspondences 
uj = (A2 + Ai) /2 and £ + 1 = a, it may appear that our N = 2 construction would coincide with 
that situation. But this is not the case. For one thing our general solution presents a periodic 
two gap structure, and is not equally spaced in its energy spectrum. Secondly, the variable G M 
takes values in the entire real line, and is not a radial variable as is r > 0. In particular, the 
integration measure over the configuration space coordinate which enters the orthonormalisation 
conditions for wave functions is thereby bound to be different as well, since indeed that example 
of a shape invariant potential applies to the three dimensional harmonic oscillator. And finally, 
the ranges of values for (^ + 1 > 1) in one case, and (—1/2 < a < +1/2) in our case, are different. 

A further difference with that example from Ref. [2] is that while the potential following 
from (j9ip is bounded below even though divergent at the finite value r = (with the range 
< r < 00), in our case among the two potential energies Vi and V2 one of these always is 
unbounded below — which of the two being dependenjfl on the sign of a or (A2 — Ai) — while 
in both cases the potential energy diverges aX u = uq. And even though all energy eigenwave 
functions remain normalisable for both potentials, these features imply specific behaviour of the 
wave functions around the singularity at n = uq, which in our case may be approached from both 
sides on the real line. 

A careful consideration of the complete solutions in ([69]) and ([70|) finds that nothing parti- 
cular occurs for the "odd" states '(/'i,2p+i('u) or V'2,2p+i('u)- Given the range (—1/2 < a < +1/2), 
these wave functions remain finite at n = uq, and in fact even vanish at that point (which, of 
course, is consistent with the sign issue for the two domains Vj^ and f_, and the fact that the 
potentials are even in (n — uq) while the "odd" states are odd functions of that variable). Given 
that the generalised Laguerre polynomial which is involved, L^^'^^°^\v), is of order p in v and 
thus of order 2p in [u — uq), these wave functions possess a total of {2p + 1) zeros (all of integer 
degree except for the one at u = uq which behaves nonanalytically as |ti — tio|^^") as expected 
of states at an odd excitation level, n = 2p + 1. However this contrasts with the states at even 
excitation level, n = 2p, which are even in the variable (n — uq) and are thus not constrained in 
their value at the point u = uq. 

Away from the value a = which corresponds to the degenerate situation of the harmonic 
oscillator which certainly is not singular at (u — uq), from (j69p and (j70p one observes the following 
situation. For the sake of the discussion, let us consider the case a > 0, knowing that the case 
a < would correspond to permuting the roles of the two systems Hi and H2 in the analysis 
hereafter. For q > 0, the potential Vi is unbounded below at u = uq, while V2 is bounded below 
but diverges to +00 at u = uq, both potentials behaving as {u — uq)~'^ in the vicinity of that 
point. Then all wave functions V'i,2p(i^) are seen to vanish as \u — uqI" at u = uq, while all wave 
functions ip2,2piu) diverge at that point as \u — uo\~'^. Yet this divergence is still just mild enough 
for these latter wave functions as well to be normalisable over M. In other words, by vanishing 
at u = Uq, the states ipi^2p{u) avoid falling "at the bottom" of the infinite negative throat-like 
potential well since then their probability density, |V'i,2p(ii)P) vanishes at the point u = uq. On 
the other hand for the states ip2,2p{'u), their probability density diverges as |n — no|~^" at u = uq 
with a power, (— 1 < —2a < 0), just mild enough for the total probability to remain finite (and 
thus normalisable to unity) and for having a nonvanishing and finite tunneling probability across 
the potential barrier in spite of a probability density which grows infinite just at u = uq, making it 
thus possible for these "even" states to have wave functions of even parity in the variable (u — uq). 
Since the generalised Laguerre polynomials Lp ^^^^"\v) involved in these wave functions are of 
order p in v, hence of order 2p in the variable {u — uq), the wave functions V'2,2p('u) indeed possess 
2p zeros (of integer degree), as expected of states at an even excitation level, n = 2p. However 

*Vi is unbounded below for a > 0, while it is V2 when q < 0. 
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for the states V'i,2p('i*) we encounter the pecuUar situation of {2p + 1) zeros, aU of integer degree 
except for the one at u = uq which behaves nonanalytically as |tt — «o|". And yet, ah these states 
are at an even excitation level, n = 2p. 

5 Conclusions 

This work addresses and solves the first two simplest cases of a general programme outlined in 
the Introduction, within the context of Supersymmetric Quantum Mechanics and its hierarchies 
of integrable quantum Hamiltonians intertwined in pairs. Heretofore in the literature in order to 
restrict the structures of such hierarchies, the requirement of "shape invariance" has been used 
for intertwined Hamiltonians. In the present work we suggest to go beyond the next level in a 
hierarchy, and rather relate Hamiltonians separated by an arbitrary number of levels. A further 
restriction even is requiring that the related Hamiltonian systems be in fact identical, differing 
only by an overall shift in their energy spectra. The added advantage of this latter requirement 
is to induce a periodic structure of finite order N in the hierarchy related to a choice of the first 
gaps in the energy spectra to be repeated then in a periodic fashion, and a choice of potential 
energies whose energy eigenfunctions would involve specific classes of orthogonal polynomials. 
Through the methods proper to Supersymmetric Quantum Mechanics, new results in the form of 
generalised Rodrigues formulae and recursion relations for such classes of orthogonal polynomials 
would also be generated. 

Besides the latter interest of a mathematical physics character, there is also a potential real 
physics interest to this approach. Indeed, by specifying the first N gaps of some energy spectrum 
for a single degree of freedom system, it should be possible to engineer through such methods a 
quantum Hamiltonian whose energy spectrum includes these N first energy levels, then repeated 
periodically, opening the way towards approximation methods amenable to numerical techniques 
as well in order to model actual quantum physical systems. And as mentioned in the Introduction, 
one last avenue to be explored further along such ideas is that of inverse scattering and integrable 
systems. 

In the present work the simple case A = 1, in which two intertwined Hamiltonians are 
identical up to the overall shift in energy, is solved as a warm-up exercise. But the detailed and 
complete solution of the N = 2 case is the actual original content of the paper, making explicit 
in that case most of the ideas outlined above. It is established that the class of orthogonal 
polynomials related to this choice is that of the generalised Laguerre polynomials. And indeed, 
new generalised Rodrigues formulae have been identified, as well as specific recursion relations. 
At the same time new classes of integrable quantum Hamiltonians with singular potentials, that 
generalise the harmonic oscillator potential, have been identified and solved completely solely 
using algebraic techniques of Supersymmetric Quantum Mechanics. Their quantum states display 
some interesting properties related to the singularity of their potential energies. 

The next step in this programme would be to manage solving the coupled Riccati recursion 
relations for N = 3, thereby opening the way in fact to the general case with arbitrary N. 

Provided this is feasible in analytic form while a numerical approach may always be developed, 
one should expect new and interesting results relating to the different directions outlined above. 
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